AA Excitation in Proton-Proton Induced 7r°7r° Production 
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Exclusive measurements of the pp —> ppTv°Tv° reaction have been performed at CELSIUS/WASA at energies 
from threshold up to T p = 1.3 GeV. Total and differential cross sections have been obtained. Here we concentrate 
on energies T v > 1 GeV, where the AA excitation becomes the leading process. No evidence is found for a 
significant ABC effect beyond that given by the conventional i-channel AA excitation. This holds also for the 
double-pionic fusion to the quasibound 2 He. The data are compared to model predictions, which are based on 
both pion and p exchange. Total and differential cross sections are at variance with these predictions and call for 
a profound modification of the p-exchange. A phenomenological modification allowing only a small p exchange 
contribution leads to a quantitative description of the data. 



Two-pion production in nucleon-nucleon col- 
lisions connects tttt dynamics with baryon and 
baryon-baryon degrees of freedom. In the spe- 
cial case that the participating nucleons fuse to a 
bound nuclear system, there is the puzzling ABC 
effect, which stands for a low-mass enhancement 
in the isoscalar 7T7t invariant mass spectrum. Very 
recent experiments on this topic are discussed in 
terms of a AA mediated isoscalar resonance in 
the baryon-baryon system as source for this pe- 



culiar ABC effect [TEE]. 

By contrast the isovector 7T7r channel in double- 
pionic fusion behaves regularly, i.e. shows no 
ABC effect and follows the expectations from con- 
ventional t-channel A A calculations 4 . Also in 
the two-pion production to unbound nuclear sys- 
tems the ABC effect was thought to be absent. 
However, a very recent inclusive measurement of 
the reaction pp — > ppX , where X stands for ejec- 
tiles not detected in the experiment, reports evi- 



dcnce for an ABC effect also in this case [5]. 

A recent isospin decomposition jB] of the to- 
tal cross sections measured in the reactions pp — > 
pp7r + TT~ , pp —¥ ppir°ir°, pp — > pnir + ir° and pp — > 
nmt + Tr + reveal these two-pion production chan- 
nels to be dominated by excitation and decay 
of resonances. In particular the TV* (1440) domi- 
nates at energies close to threshold and the AA 
system and possibly the A(1600) at higher in- 
cident energies. The latter is supposed to con- 
tribute primarily to the pp — > nmr + Tr + and also 
to the pp — > pmr + Tr° channel. 

In view of the challenging interpretation [5] of- 
fered for the ABC effect in isoscalar irir chan- 
nels in case of double-pionic fusion and the re- 
ported evidence in the inclusive pp — > ppX reac- 
tion it appears mandatory to study the isoscalar 
7T7T production with exclusive and kinematically 
complete measurements in the case, where the 
two participating nucleons do not fuse into a final 
nuclear bound system. Among the two possible 
choices, the ppir + ir~ or the ppir°Tr° channel, the 
latter one is especially appealing, since it contains 
no 7T7T isovector contributions, only isoscalar and 
isotensor parts with the isoscalar part being the 
by far dominating one [6 

From previous work it is known that the pp — > 
ppir + ir~ and pp — > ppir°ir° reactions in the near- 
threshold region are well understood as being 
dominated by excitation and decay of the Roper 
resonance |7I8I9I10| . At higher energies theoreti- 
cal calculations [7] predict the t-channel AA ex- 
citation to play the dominant part. These cal- 
culations are compared in Figs. 1-5 with the 
differential and total cross section data for the 
pp — > ppir°ir° reaction obtained in this work. 

Since there have been no exclusive measure- 
ments of the pp — > ppir°ir° channel in the en- 
ergy region of interest, we have carried out a sys- 
tematic program of exclusive two-pion produc- 
tion measurements in pp collisions from thresh- 
old up to T p = 1.36 GeV using the WAS A detec- 
tor [11] with the hydrogen pellet target system at 
the CELSIUS storage ring of the The Svedberg 
Laboratory in Uppsala. The detector has nearly 
full angular coverage for the detection of charged 
particles and photons. The forward detector con- 
sists of a thin-walled window plastic scintillator 



hodoscope at the exit of the scattering chamber, 
followed by straw tracker, plastic scintillator quirl 
and range hodoscopes, whereas the central de- 
tector comprises an electromagnetic calorimeter 
consisting of 1012 Csl (Na) crystals, and in its in- 
ner part a plastic scintillator barrel surrounding 
a thin-walled superconducting magnet containing 
a mini drift chamber for tracking. 

Neutral pions are reconstructed from photons 
detected and identified in the central detector. 
Protons are detected in the forward detector and 
identified by the AE-E technique. Since the for- 
ward detector cone does not cover the full kine- 
matic angular range for protons at high incident 
energies, the detection efficiency for protons at 
medium center-of-mass (cms) angles is reduced. 
This introduces systematic uncertainties in par- 
ticular in the proton angular distribution. Since 
due to the identity of the two incident particles, 
the angular distributions have to be symmetric 
about 90°, the observed asymmetries about 90° 
(see Figs. 3 and 4) may hence serve as a mea- 
sure of such systematic errors in the data. We 
estimate the systematic uncertainties due to this 
deficiency of full phase space coverage by using 
various model calculations in the Monte Carlo 
simulations of the detector response and accep- 
tance corrections. The estimated systematic un- 
certainties are shown by dark-shaded histograms 
in Figs. 2-5. 

The absolute normalization of the data has 
been achieved by a simultaneous measurement of 
clastic scattering and/or single pion production, 
for which the cross sections are known. Since in 
particular the single-pion production cross sec- 
tions are not known better than to an accuracy of 
20 %, this uncertainty transmits also to the cross 
sections deduced for the pp —¥ ppir°ir reaction. 
For details of the data analysis see, e.g. Refs. 

m- 

Experimental results for the low-energy range 
T p < 1 GeV have been published [TU] in connec- 
tion with the discussion of the properties of the 
Roper resonance. In addition close-to-threshold 
results from previous PROMICE/WASA mea- 
surements are given in Ref. [12]. 

The total cross section data from this work 
have been published already [6] in connection 
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Figure 1. Total cross sections for the pp — > 
ppir°Tr° reaction. Data are represented by open 
(bubble chamber data, Refs. }13|14j ) and filled 
symbols (WASA data, Refs. [EG]). The dot- 
ted lines show the original calculation of Ref . [7J . 
The dashed, dash-dotted and solid lines are calcu- 
lations with modifications described in the text. 



with the isospin decomposition of two-pion pro- 
duction data. They are shown in Fig. 1 to- 
gether with previous bubble-chamber results [131 
[14] . The total cross section keeps rising from 
threshold up to T p « 1 GeV, where it levels off 
and proceeds only slowly rising until 1.2 GeV. 
Thereafter it continues steeply rising until 1.5 
GeV, where it finally levels off again - see Fig. 
3 in Ref. [6 . As has been demonstrated [6], the 
low-energy structure is due to the Roper reso- 
nance, whereas the renewed rise at higher ener- 
gies can be associated with the dominance of the 
AA excitation. 

In Figs. 2 - 5 we exhibit a selection of eight 
distributions, which are most significant with re- 
gard to their physics content. Note that for a 
four-body reaction with unpolarized beam and 
target there are seven independent single differ- 
ential distributions. For T p — 1.0, 1.1, 1.2 and 
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Figure 2. Distribution of the 7r°7r° invariant mass 
M^o^o (left) and the 7r°7r° opening angle 6 v o v o 
(right) for the pp — > ppTr°ir° reaction at beam 
energies T p = 1.0, 1.1, 1.2 and 1.3 GeV (from 
top to bottom). Solid dots represent the ex- 
perimental results of this work. The light-shaded 
areas denote phase space distributions and dark- 
shaded histograms systematic uncertainties. The 
dotted lines show the original calculation of Ref. 
[7j. The dashed, dash-dotted and solid lines are 
calculations with modifications described in the 
text. All calculations are normalized in area to 
the data. 
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Figure 3. Same as Fig. 2 but for the distributions 
of the invariant masses M p7I .o (left) and the M pp7T a 
(right) for the pp — > ppir°Tr° reaction at beam 
energies T p = 1.0, 1.1, 1.2 and 1.3 GeV (from 
top to bottom). 



Figure 4. Same as Fig. 2 but for the distri- 
butions of the cms angles 9 P (left) and the 0„.o 
(right) for the pp — ► ppir a ir° reaction at beam en- 
ergies T p = 1.0, 1.1, 1.2 and 1.3 GeV (from top 
to bottom). 
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Figure 5. Same as Fig. 2 but for the distributions 
of the angle of the tt°tt° system in the cms S^o^o 
(left) and the tt° angle in the ir°ir° subsystem 
G^o 71 " (right) for the pp — > ppir°ir° reaction at 
beam energies T p — 1.0, 1.1, 1.2 and 1.3 GeV 
(from top to bottom). 



1.3 GeV differential distributions are shown for 
the invariant masses M v o w o, M p7V o, M ppw o , the 
opening angle between the two pions S-^o^o, the 
proton angle G p , the ir° angle G„.o, the angle of 
the 7r 7r° system G^o^o - all in the cms - as well 
as the 7T° angle in the 7r°7r° subsystem (Jackson 
frame) . 

In Figs. 1-5 the data are compared to pure 
phase space distributions (light-shaded areas in 
Figs. 2 - 7) as well as to calculations of Ref. 
[7] with and without modifications as will be dis- 
cussed in the following. As a convention we show 
all theoretical model distributions normalized to 
the experimental total cross section, i.e. to the 
same area in the differential distributions. This is 
because we are interested here in the shape of the 
differential distributions. The comparison with 
the absolute total cross sections is done in Fig. 1. 

We see that many of the experimental differ- 
ential distributions are not far from the phase 
space distributions. With regard to angular dis- 
tributions significant deviations from isotropy are 
observed only for the G p and ^o„o distributions. 
The first one is largely characterized by the t- 
channel exchange as demonstrated in Ref. [5], 
whereas the latter one is strongly correlated with 
the M^o^o spectrum as discussed in some detail in 
Ref. [9 . In fact, we observe some deviations from 
phase space also for the M n o„o spectrum as will 
be discussed below. However, really large devi- 
ations from phase space are observed in spectra, 
which are correlated with A excitation: M p7r o, 
M pw o n o and M pp7r o. The peaks in all these in- 
variant mass spectra build up increasingly with 
increasing energy and reflect the increasing ex- 
citation of the AA system, the resonance pole of 
which is reached at T p « 1.3 GeV, i.e. the highest 
energy considered in this work. The M pp spec- 
trum, a sample of which is shown in part in Fig. 
6, is kinematically complementary to the M^o^o 
spectrum. 

Next we confront the data with model predic- 
tions [7 of the Valencia group and subsequent 
modifications of the original calculations. The 
dotted lines in Figs. 1-5 show the original pre- 
dictions, which though renormalized in area give 
huge deviations from the measured distributions - 
in particular in the M n o w o and S w o n o spectra. For 



the M^o^o spectrum these calculations predict a 
kind of double-hump structure with a large en- 
hancement at high invariant masses, which is ab- 
sent in the data. With regard to the 7r°7r opening 
angle S^o^o these calculations predict a preferen- 
tial antiparallel emission of the two pions, which 
again is not supported by the data. 

The predicted double-hump structure is remi- 
niscent of similar predictions for the ABC effect 
in the doublc-pionic fusion. There a double-hump 
structure has been predicted [15116117] based on a 
^-channel A A excitation — with some of such cal- 
culations |18U9j emphasizing a particular strong 
high-mass enhancement. By contrast exclusive 
and kinematically complete measurements of the 
isoscalar double-pionic fusion find only a huge 
low-mass enhancement, but no or no significant 
high-mass enhancement [T|2 20 . As already men- 
tioned above our data for the M x o T o spectrum are 
in qualitative agreement with phase space, i.e. we 
observe neither a spectacular low-mass enhance- 
ment (ABC effect) nor a spectacular high-mass 
enhancement. 

In order to shed some light into the failure of 
the theoretical predictions, we readjust these cal- 
culations step by step first by implementing the 
knowledge accumulated from the studies of the 
near-threshold region (T p < 1 GeV), which is 
governed by excitation and decay of the Roper 
resonance. From the analysis of the data at T p 
= 0.775 and 0.895 GeV [TO] we obtain a value 
for the relative branching between the decay via 
the A resonance, i.e. N* — > Air — > N(tttt)i = i = o 
and the direct decay N* — » N(Tnr)i = i =0 , which is 
four times smaller than that quoted in PDG [21] 
and used in Ref. [7J, but in agreement with a re- 
cent analysis of data on pion- and photo-induced 
pion production on the nucleon [22]. Updating 
the model calculations with this new branching 
for the Roper decay (dashed lines in Figs. 1-5) 
leads not only to a quantitative description of the 
data for T p < 1 GeV [TO], but also to a consid- 
erable improvement in the description for T p > 1 
GeV, though there are still substantial deficien- 
cies in the M^o^o and S^o^o distributions. 

Next we readjust the total strength of the 
Roper excitation. We know from the isospin de- 
composition 6 of the total two-pion production 



cross sections that for T p > 1 GeV the excita- 
tion of the Roper resonance comes out much too 
strong in the Valencia model calculations: Hence 
we readjust the strength of the Roper excita- 
tion according to the isospin decomposition re- 
sult. That way we force the calculations to also 
reproduce the total cross sections. The outcome 
of this modified calculation is shown in Figs. 1 - 
5 by the dash-dotted lines. We now get a good 
description of the data at T p = 1.0 GeV. How- 
ever, deviations from the data still increase with 
increasing beam energy. In particular, because 
the last modification increased strongly the dom- 
inance of AA excitation for T p > 1 GeV, we 
obtain pronounced double-hump and parabolic 
structures, respectively, in the M^o^o and S^o^o 
distributions. 

Since these failures are now intimately con- 
nected to the treatment of the AA excitation, 
further improvements have to be sought in a mod- 
ification of this process. As demonstrated already 
in Ref. [9], double-hump structures in M w o w o 
spectra are generated by a k\ ■ k^ term, where 
fciand ki denote the 3- momenta of the two emit- 
ted pions. In the description of the A A process 
we find indeed such a k\ ■ k^ term associated with 
the p exchange, see equation (A. 10) in Ref. [7J. A 
closer inspection of eq. (A. 10) shows the follow- 
ing problem: the p exchange provides isotropic 
angular distributions for pions as required by the 
data, however, provides also a k\ ■ ki term at 
variance with the data. On the contrary, the 
pion exchange contains no k\ ■ ki term, but in- 
volves a strong angular dependence for the pi- 
ons - in disagreement with the data. Whereas 
pion-exchange appears to be straightforward the- 
oretically and well established in the description 
of pion-production processes, stringent tests are 
missing for the p-exchange. Also since the p- 
cxchange part in the i-channel AA description 
of Ref. [7] is the by far dominating part, we have 
to seek for a solution in the ansatz for the p ex- 
change. 

In a recent theoretical work on two-pion pro- 
duction Cao, Zou and Xu [23 find that, in con- 
trast to Ref. [7] , p-exchange is small compared to 
pion-exchange. Triggered by this finding and in 



view of our problems to describe the data within 
the context of the ansatz equation (A. 10) in 
Ref. [7], we investigated the possibility to achieve 
a reasonable description of the data by varying 
the coupling strength of p exchange. For simpli- 
fication, compared to eq. (22) of Ref. [7], we use 
the longitudinal and transversal exchange contri- 
butions to be given purely by n- and p-exchange, 
respectively. Also for sake of simplicity we drop 
the monopole meson form factors. 

In order to have the description relativistically 
more appropriate we take the momenta of the 
pions from the A decay in the corresponding A 
systems and correct them by the Blatt-Weisskopf 
barrier factors as given by Pilkuhn ^3] and used 
also in Refs. [15123] . These relativistic correc- 
tions are, of course, negligible close to thresh- 
old, but lead to sizable and significant changes 
at higher energies. In particular the calculated 
total cross sections no longer keep steeply rising 
beyond the resonance pole as in Ref. [7J, but go 
into saturation as shown in Ref. [23]. The latter 
behavior agrees with the trend of data at high 
energies. 

In fact, we find a quite reasonable description 
of the data by reducing the coupling strength of 
the p exchange (A. 10) of Ref .^7 by an order of 
magnitude. The solid lines in Figs. 1-5 show our 
calculations with the p coupling strength reduced 
by a factor of 12 and reversed in sign. Note that 
now the p-exchange is only a small correction to 
the leading 7r-exchange, but this correction leads 
still to a sizable improvement in the description 
of the differential data. The description of the 
data is still far from being perfect, but at least 
all major features of the data are reproduced. 

Fig. 1 shows the energy dependence of the to- 
tal cross section of the pp — > ppir°ir° reaction. As 
mentioned in the introduction the striking feature 
is the slow rise of the cross section between T p = 
1.0 and 1.2 GeV, which is reproduced neither by 
the calculations of Ref. [7] nor by the more re- 
cent ones of Ref. [23] . The primary reason for that 
failure is that in both calculations the Roper ex- 
citation keeps rising beyond T p = 1.0 GeV, which 
is at variance with the finding from the isospin 
decomposition [BJ. As pointed out above the dif- 
ferential cross sections require the Roper excita- 



tion to be cut down, too, in accordance with the 
isospin decomposition result. Hence, accounting 
for the latter and using the modified AA descrip- 
tion we succeed in obtaining a description for the 
total cross section, which is in quantitative agree- 
ment with the data. 

The comparison to the data shows that the p 
exchange contribution as treated in Ref. |7|18) 
needs substantial modification. This is true both 
for the two-pion production to unbound nuclear 
systems and for the double-pionic fusion. No sig- 
nificant high-mass enhancement is observed in ei- 
ther case, which would be a signature of a domi- 
nant p-exchangc. 

In the double-pionic fusion experiments a huge 
low-mass enhancement (ABC-effect) is found in- 
stead, which has been associated with the for- 
mation of an isoscalar resonance via a AA door- 
way [? . Such an isoscalar resonance can not con- 
tribute to the isovector pp — > ppir°ir° channel dis- 
cussed here. The small low-mass enhancement 
visible in the data can be fully associated with 
the conventional t-channcl AA excitation. 

In Ref. [5] the appealing idea has been put for- 
ward that though there is no bound state in the 
isovector NN system one could look for a qua- 
sibound 2 He double-pionic fusion process in the 
reaction pp — > ppinr by requiring that the emitted 
protons have very small kinetic relative energies. 
In Ref. 5 , which presents and discusses COSY- 
ANKE data for the inclusive reaction pp — > ppX, 
this requirement was achieved by the condition 
M pp < 2m p + 3 MeV. These data — covering 
just the forward angle region — exhibit a pro- 
nounced low-mass enhancement in the pp missing 
mass spectrum, which is equivalent to the asso- 
ciated 7T7T invariant spectrum. In our experiment 
we cover practically the full phase space of the 
two-pion production reactions. Though we have 
accumulated quite some statistics of events dis- 
tributed over the full phase space, we lose nearly 
all events, if we apply the above M pp cut. There- 
fore, in order to have acceptable statistics, we re- 
lax the above constraint to M pp < 2m p + 10 MeV. 
In Fig. 6 we show for T p = 1.3 GeV — the energy, 
where the AA process is most obvious in our data 
— the Mpp spectrum together with the indicated 
cut as well as the M^o^o, M p7V o and 6^(1^0 spectra 
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Figure 6. Differential distributions at T p = 
1.3 GeV for the 2 He scenario. Top left: M pp 
spectrum for the full phase space (only low-mass 
part shown). The distributions for M^o^o (top 
right), M p7r o (bottom left) and O^o^o (bot- 
tom right) are plotted with the 2 He condition 
Mpp < 2m p + 10 MeV (vertical line in the M pp 
spectrum). The drawn lines represent calcula- 
tions as described in the caption of Fig. 2 and 
text. The shaded areas show the phase space dis- 
tributions. 



resulting from this cut. In the M pp spectrum we 
observe an indication of the low-mass enhance- 
ment produced by the pp final-state interaction, 
which is also included in the theoretical calcu- 
lations. In order to focus on small masses we 
plot this spectrum only in the low-mass range in 
Fig. 6. We see that the cut M pp < 2m p + 10 
MeV is still within the region dominated by the 
pp final-state interaction ensuring thus relative s- 
waves between the two protons. Within the lim- 
ited statistics the data in the constrained M n o„o 
spectrum exhibit some high-mass enhancement 
as expected from the pioneering work of Risser 
and Shuster [15] concerning the production of an 



isoscalar pion pair via the AA excitation in the 
double-pionic fusion. The constrained M p7r o spec- 
trum exhibits the A peak ensuring the AA pro- 
cess to be the dominating process in the 2 He case, 
too. 

The solid lines, which present the modified the- 
oretical description of the pp —¥ ppir°ir° reaction, 
give a reasonable description for this constrained 
phase space scenario. On the contrary, the bro- 
ken lines exhibit a very pronounced double-hump 
structure in the M„o„o spectrum due to the dom- 
inance of the p-exchange in these calculations. 
This behavior is clearly at variance with the data. 
On basis of the model calculation we also ensure 
that the change in the cut from a 3 MeV to a 
10 MeV range does not change the results qual- 
itatively. The 8£™ angular distribution points 
to some sizable d-wave contribution, which is 
in support of the respective ansatz in Ref. [5]. 
The increased d-wave contribution as compared 
to the results for the full pp — > ppTr°ir° reaction 
is not unexpected, since the 2 He cut forces the 
two emerging protons to be in relative s-wave. 
The decay of the two excited A states results in 
a double p-wave emission of the two pions. This 
in turn leads to s- and d-waves in the emitted 
ppir°ir° system. 

Finally we compare our calculations to the 
COSY-ANKE data in Fig. 7 by restricting the 
calculations further to the angular range covered 
by the ANKE experiment, which is cosQ pp > 0.95 
with Q pp = 180° - 6,^0 and 9 l ° b < 12°. We 
again obtain an essentially quantitative agree- 
ment with the data at T p = 1.1 and 1.4 GeV us- 
ing the modified model description, whereas the 
p-exchange dominated calculations again give a 
double-hump structure in vast disagreement with 
the data. 

Summarizing, we have presented first exclusive 
and kinematically complete measurements of the 
AA system excited in the pp — > ppir°ir° reaction. 
The data are well described by a conventional t- 
channel AA calculation, where the p exchange 
contribution is strongly reduced compared to that 
in Ref. |7]. The same holds for the subset of 
data, which corresponds to the double-pionic fu- 
sion to a quasibound 2 He nucleus. No evidence 
is found for a significant low-mass enhancement 
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Figure 7. M^o^o spectra at T p — 
for the 2 He condition M pp < 2m p 
ANKE angular constraints cosQ pp > 0.95 with 
O pp = 180° - 9^0 and Q l p ah < 12°. The solid 
dots represent the ANKE data [5] and the drawn 
lines calculations as described in the caption of 
Fig. 2 and text. The shaded areas show the phase 
space distributions. 



in the M pi o^o spectra (ABC effect) beyond that 
given by the conventional t-channel AA excita- 
tion process. The new data sets should serve as a 
significant test case for theoretical treatments of 
multipion production. 
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